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Abstract 

An analytical approach based on the parametric representation of the 
wave propagation in nonuniform media was considered. In addition to 
the previously developed theory of parametric antiresonance describing 
the field attenuation in stop bands, in the present paper the behaviour of 
the Bloch wave in a transmission band was investigated. A wide class of 
exact solutions was found and the correspondence to the quasi-periodic 
Floquet solutions was shown. 
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Introduction 

In this paper we concentrate on the analytical theory of linear propagation in 
nonuniform media based. The difficulty of this problem comes out even in the 
one-dimensional (ID) the ordinary differential equation (ODE) 

w"(x) + q 2 (x)w(x) = 0, q(x) = —n(x), (1) 
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has no explicit solution except for some rare model functions q(x). The above 
ODE describes ID harmonic waves ~ exp (— iwt) propagating in a nonuniform 
dielectric medium with gradually varying refraction index n(x); c is the speed 
of light in vacuum and / denotes the differentiation with respect to x. 

Real-valued solutions of ODE (fTJ) correspond to standing waves carrying 
no energy while complex-valued ones describe travelling waves. For q = const 
both types are "sinusoidal" (proportional to sin/ cos qx or exp (±iqx)). Evi- 
dently, variations of the parameter q(x) may change their amplitude and wave 
form. Moreover, for a periodic modulation of the refraction index n(x) the 
parametric resonance/antiresonance (i.e., the existence of exponentially increas- 
ing/decreasing solutions) may occur. 

According to the Floquet theorem [TJ [2] , for any general periodic modulation 
q (x + x) = q( x ) ( or equivalently n (x + x) — n ( x ))t where \ is the period, there 
exists the following solution of the equation (TTJ) (Bloch wave) : 

w(x) = w(x) exp (ifix) , (2) 

where w(x) is a periodic function and the characteristic exponent /x can be ei- 
ther (i) real or (ii) purely imaginary [2], The former case (with the minus sign) 
corresponds to a parametric antiresonance (wave attenuation) in the stop bands 
("forbidden zones") of the periodic structure (multilayer interference mirrors, 
Bragg optical waveguides) ; in other words it provides localization of the prop- 
agating mode by the periodic dielectric cladding. The latter case describes a 
periodic modulation of the travelling waves. 

In many applications the inverse problem (a kind of optimal control) arises 
when one needs to find such a dependence n(x) that provides the maximal 
value of the characteristic exponent \i. In our previous papers O HJ [5] a sim- 
ple analytical theory of parametric (anti)resonance was developed and explicit 
formulae for field attenuation in stop bands was derived. However, for the sake 
of completeness, we need also to know the behaviour of the Bloch wave in a 
transmission band where the characteristic exponent is purely imaginary. The 
first attempt to apply the method of phase parameter to the transmission bands 
has been reported in [B]. 

In the present paper, being a rectification of [BJ, we develop a general 
approach based on the parametric representation of the wave propagation in 
nonuniform media (both in transmission and stop bands). A physically correct 
definition of phase is put forward for a general case of linear wave propagation in 
a nonuniform medium. A wide class of exact solutions to the ID wave equation 
(fT) is found and the correspondence to the quasi-periodic Floquet solutions is 
shown. 

1 Phase parameter for the real case 

Let us recall the results obtained in the previous papers of the first author (AP) 
concerning the parametric representation of the fundamental system of solutions 
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to the wave equation (JTJ) . If we introduce an admittance function 

w'(x) 

2/W = TTTT- 
q(x)w(x) 

i.e., the solution of the equation ((I) can be represented as follows: 



w(x) = wq exp 



y(x)q(x)dx 



(3) 



(4) 



where wo is an integration constant, then it is easy to observe that the primary 
wave equation ([T]) can be equivalently rewritten in the following form: 



q{x)y'{x) + q'{x)y{x) + q 2 (x) [l + y 2 (x)] = 0, 



i.e., 



y'(x)dx 



y{x)q' (x)dx 



x - x, 



q{x)[l + y 2 (x)} J q 2 (x)[l + y 2 {x)] 

where xo is an integration constant. 

If q{x) = qo is a constant, then the above equation (j6)) reads 

dy 



(5) 
(6) 



1 



xq — x. 



qo J 1 + y 2 

We see that the integral in ([7]) can be easily integrated with the substitution 

dip o I 



(7) 



y = ctg-0, dy 



• 2 7 ' 

sin ip 



sin 2 ip ' 



which gives us the solution 

ip = ctg~ 1 y = ip + q (x- x ) = q (x- x ) , 

where ip® is an integration constant and xq — xq — ipo/lo- 

Therefore, the solution (U) can be rewritten in the equivalent form: 

w(x) = w exp | g J ctg [qo (x - x Q )] dx 
dsin [q (x - x )] 



(8) 
(9) 



wq exp 



sin 



[q (x - x )] 



— w sin [q (x - x )] , (10) 



where wq is new integration constant. 

More generally, for an arbitrary function q(x) we can, by analogy to the 
above reasoning, define a new variable ip{x) (phase parameter) as follows: 



Then the equation ([5]) reads 

_ ff(x)-0'(x) 

sin 2 ip(x) 



y(x) = ctgtp(x). 



q'(x) ctg ip(x) 



sin 2 ip(x) 



(11) 



(12) 
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i.e., 



2q(x) 



(13) 



If there exists the inversion x = X(ijj), then we can write w(x), y{x), and 
q(x) as functions of ip, i.e., 



w [X(i/j)] = y [X{4>)] = Y(tJj) = ctg V, q [X(i/>)] = Q(# 

Then the equations Q and ([6]) can be rewritten as follows: 



X(^) 



wo sin ip exp 



cos ip dtp 



x 



Q«0 . 9 / 

sm 2^y of^y. 



Q 2 W>) 



(14) 



(15) 



(16) 



Here and below dots denote the differentiation with respect to ip. For any 
arbitrary chosen function Q(ip) we obtain explicit expressions for W(i/j) and 
X(ip), therefore, for w(x) = W [ip{x)\. 

If we make a substitution Q(ip) = qoexpG(ip), where qo is some constant, 
then the above equations ([15]) and (fT6|) will take the following form: 



W{ip) 



xty) 



wo sin ip exp 



xo H 

9o 



G(ip) cos 2 d?/' 



cty> 1 f G(ip) sin 2^ dip 

expG(^) ~ 2 7 expG(V>) 



(17) 



(18) 



In such a way we construct a closed-form parametric solution which reveals 
some quantitative relations inherent in wave propagation in nonuniform media, 
not obvious in classical solutions of the corresponding differential equations pQ . 
In particular, for any periodic refraction index n(x) = n(x+x) defined implicitly 
by a Fourier series, i.e., 



G(tp) = a + ^ (<i2m cos 2mip + 6 2m sin 2mip) , 



(19) 



where ao, 0b2mi ^2m are some constants, we obtain a Floquet solution ^ with the 
minus sign in the exponent and the characteristic exponent /i = v/x with the 
explicit formulae for the period A and the attenuation per period v [3J 01 [B] : 



X 



9o Jo 



exp [—G(ip)] sin 2 ipdip, v 



G(ip)sin2ipdip. (20) 



These analytical relations, giving the very simple description of the wave field 
attenuation in a periodic structure, are useful for the optimal design of multi- 
layer mirrors and Bragg fibre claddings [3J. However, from the theoretical point 
of view this solution remains incomplete until a similar parametric representa- 
tion is found for the propagating waves in the corresponding transmission bands 
of a periodic medium. 
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2 Complex case 



For a complex wave also is possible to define a phase parameter ip(x), which 
obviously must be a homogeneous function of w(x) and w'(x). 
Let us observe that using the identity 



y(x)+i ctg ip + i cos ip -\-ishxip 
y(x) — i ctg ip — i cos i\) — i sin ijj 

Eq. (If f I) can be equivalently rewritten as follows: 



exp (2iip) 



tp(x) = ctg 1 y(x) 



I y(x) + i 



f u/(:r) + iq(x)w(x) 
In ■ 



2i y{x) — i 2i w' (x) — iq(x)w(x) 



(21) 



(22) 



Then let us define the quasi-phase parameter ip(x) of a complex wave func- 
tion w(x) as follows: 



1 w/(x) +ig(cc)w(x) 

= ttt In -TTT\ ■ < \ *( \ 

li w (x) — iq\x)w (x) 



q(x) + 



q'(x) Re[y(x)} 
q(x) \y{ X ) + if 



dx. (23) 



The second part of Eq. (|23j) can be obtained from the first one by the direct 
calculation of the integral representation of the logarithm, i.e., 



Re 



d{w' + iqw) 
w' + iqw 



= Re 


[t/ 







f i(qw' + q'w) 
w' + iqw 



dx 



(24) 



and now, using Eq. ([T} and the facts that Re (iz) = — Im(z), Im(iz) = Re(z), 
we finally obtain that Eq. (|2"41 can be rewritten as follows: 



Im 



iq (w' + iqw) + iq'w 



dx 



q+ —Re 



y(x) + i 



w' + iqw 

The complex-valued admittance y(x) as a function of tp reads 



dx. 



y [Xty)] = Yty) = 
It follows directly from Eq. (|23p that 



1 



X(<P)Q(iP)W(tP) 



G{1>) Re[Y(ij)} 



dx Xty) ' X{i!) \Y(^)+i\ 2 ' 

then Eq. (1) can be rewritten as a pair of nonlinear differential equations 

G Re Y \ 



X 



Y 



1 



y\ \Y + i\ A ) 
G Im Y [i (Y 2 -l) - 2Y] 
\Y + i\ 2 



l + Y 2 



(25) 
(26) 
(27) 

(28) 
(29) 
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where we use the fact that Q(tp) = G(ip)Q(ip). The second equation above is 
obtained if we insert w"(x) = [q(x)w(x)h(x)]' into Eq. ([1]), i.e., 



w" [X($)\ = Q 2 W 



Y + YG 
QX 



Y 2 



-Q 2 W = -q 2 w. 



This provides us also with the compatibility condition (cf. Eq. (J5J) 

Y + YG+(1 + Y 2 )QX = 



(30) 



(31) 



imposing constraints on choosing the complex admittance Y(ip). 

It is interesting to note that according to Eqs. (|2"5|) and pip we can write 
that 



• Y + GY GReY 
— : TTo~~ — 1 — 



l + Y 2 



Y + i 



or equivalently 



Y = G 



(1 + Y 2 ) Re Y 



\Y + t\ 2 



Y 



(l + Y 2 ) 



This is another way to write Eq. 



Let us also note that 



\y(x)+if = \y(x)\ 2 + 2Im[y(x)]+l. 



(32) 



(33) 



(34) 



If the functions Q(ip) and Y(ip) are given, then using Eqs. (US]), O, P5|l. 
and (|31[) the following expressions for X(ifj) and the complex- valued wave func- 
tion W(ip) can be written: 



W(ip) 



GReY\ dtp 



\Y + i\ 



Q ' 



exp 



QXYdtp 



u>o exp 



(35) 
(36) 



3 7Z, y~ variables 

For the sake of convenience, let us denote Y = lZexp(iy) and separate real and 
imaginary parts of Eq. (|29| (or equivalently Eq. (1331) ). then as a result we obtain 
the following pair of nonlinear differential equations: 



n = -GS{K,y)[21l + (l + K 2 )smy] ~ (l + K 2 )cosy 



y = 



-GR + (K 2 + 1 
K 2 -l 



GC(K,y) - 1 



cosy, 



n 



GC(K,y)-l 



siny, 



(37) 
(38) 



G 



where 

s(K,y) = TTwri v > C ( K ^)- 1+ Jr^ v (39) 

1 + 7c 2 + 27c sin y 1 + 72 ^ + 272 sin J 7 

In new variables we have that 

Re F(V>) = 72(V>) cos^(^), Im Y(tp) = K(ip) sin^(V'), (40) 
and Eq. (|34[) takes the following form: 

+ i| 2 = l+ft 2 (^)+27t(V0sin;y(# (41) 
Then the functions S (7c, y) and C (72., y) can be defined as follows: 

\Y[ip)+t\ \Y(tp)+t\ 
with the compatibility condition 

S 2 (72, y) + C 2 (72, y) = '^y 1 ' . (43) 

Then Eq. ([29]) reads that 

(72 + i72y) = GS (72, y) [t (72 2 e 2 ^ - l) - 272e^] - (l + 72 2 e 2l3; ) . (44) 
The real and imaginary parts of Eq. (|44]) can be written as follows: 



72cosy- 723>siny = -GS (72,y) (1 + 72siny) 272cosy 

- (l + 72 2 cos2y) , (45) 

72siny + 723>cosy = -GS (72, y) (l + 272 siny - 72 2 cos 2y) 

-72 2 sin2y. (46) 



Then multiplying Eq. (|45|) by cos y and Eq. (|46|) by sin y and taking their sum 
we obtain Eq. (|37|) and multiplying Eq. (j45]) by sin y and Eq. (|46|) by cos y and 
taking their difference we obtain Eq. (f38j) . Additionally excluding from Eqs. 
(|4"5"|) and (|46p the term with the derivative G we obtain that 



-r^(v v\ = ncos y- n y sin y + 1 + n2 c ° s 2 ^ 

1 272 cos y (1 + 72 sin y) 

72 sin y + TZy cos y + 72 2 sin 2y 
1 + 272 siny- 72 2 cos 2y ' 

i.e., equivalently 



(47) 



72 cosy + (1 + 72 2 ) + 272siny 



(1 - 72 2 ) = 723> [272 + (l + 72 2 ) siny] . (48) 
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This is the compatibility condition which is satisfied by TZ(i/j) and y(ip). 

Let us note that Eqs. ([37]) and (|38)l can be equivalently rewritten as follows: 



K + GK 
K 2 + l 

ny 

Tl 2 - 1 



GC(n,y)-i 

GC(K,y)-l 



cos J 7 , 
sin^. 



(49) 
(50) 



In a general complex case Eqs. ([57| and (155)) can be integrated with respect 
to yty): 



ywo 

Really, let us denote 



arcsm < — _ , , , ex P 



J(n,y) 



-2 



i + n 2 {^) 



ft sin ^ 



(51) 



(52) 



ft 2 + 1 

Then using Eqs. ([57)) and (I3"8l we can calculate its derivative with respect to ip: 

x (i - n 2 ) 11 sin y + (i + TZ 2 )ny cosy n 2 -i- , , , s 
^ (W y) = ^ \, 2 — — — - = ^t-^gj (n, y) . (53) 



{n 2 + \y 

We see that Eq. ([55)) can be easily integrated: 

j(i>) = Jin(<<p),y(i>)} = £(Tp) 

where 



TZ 2 + 1 



£(ip) = exp 



WWTi G( W 

and we used the fact that 



= Q(ip) exp 



-2 



Finally, we obtain that 



and this leads us directly to Eq. ([51 



s in yw = i±#^)^) 



(54) 
(55) 

(56) 
(57) 



4 Second-order nonlinear differential equation 

The functions S(ip) = S [TZ(tp) , y (tp)] and £(ip) (or l 7(V') equivalently) are con- 
nected by the following relation: 



s(V0 



1 + 2£(V>) 



(58) 
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We can also calculate the following derivatives with respect to ip: 



Kcosy 



s{n,y) 



n 2 + i 



ft cosy 



(1 - ft 2 ) Ttcosy - (i + n 2 ) nysmy 



n 2 + 1 



ft 2 -l 
K 2 + 1 

c(n,y) 
n 2 -i 



(k 2 + ly 
\ + 2GS{n 1 y) 

ft sin y 



K cosy 



K 2 + l 



n 2 + i 



n 2 



-is(n,y) + G 



s(n,y) 



Kcosy 



n 2 + i 



(59) 



(60) 



For the function C(i[>) — C [TZ(ip),y(ip)} we obtain a nice nonlinear second- 
order differential equation 



C(V0 + 4C(V) = 



C 2 (^)+4C 2 W-1 



(61) 



with the eigenfrequency 2 and modulation determined by the variable refraction 
index n{x) = uq exp [G \^(x)]\, where uq — (c/w) qo- 

Really, using Eqs. (|57|) and we can calculate the first and second deriva- 
tives of the function C(ip) as follows: 



C 



(i-ft 2 ) n cosy - [(i + n 2 )smy + 2ii\ ny 



(1+ft 2 + 2ft sin XT 



ft 2 -l 



l+ft 2 + 2ftsin;y 



and 



C 



-25-2 



1-25 



1 



1 + K 2 +2ftsin^ 



l+ft 2 + 2ftsi n ;y 



(62) 



(63) 



Calculating the auxiliary derivatives 

(1 - K 2 ) ft sin ^ + (1 + ft 2 ) ft3>cos^ 



S 



(l+ft 2 + 2ftsin^) 2 



= GS- 2GS 2 - 2GS 



1 



l+ft 2 + 2ftsiny 



and 



1 



l + ft 2 + 2ftsiny 



= 2C + 2GS 2 + 2GS 



l+ft 2 + 2ftsin^ 



(64) 



(65) 



we see that their second and third terms are eliminating each other and then 

C = -AC - 2GS. (66) 
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Let us also note that 

^ 2 2 _ (ft 2 - 1) 2 + 4ft 2 cos 2 y _ (ft 2 + 1) 2 - 4ft 2 sin 2 y 
+ ~ (ft 2 + if (1 + 2jf ~ (ft 2 + l) 2 (1 + 2J) 2 

= 1 " 4J2 2 =^ = 1-4S- (67) 
(1 + 2.7) 2 1 + 2J V ; 

Finally, excluding 5 from Eqs. (j6"6")l and (1571) we obtain Eq. (jSTj) . 
There are two ways of constructing sought parametric solutions: 

(i) to define G(ip) and then solve Eq. (fUTj) with respect to C(V>) or 

(ii) to define C(V0 and then find G(tp) by integration: 

G W = 2f. d M +4C M #. (68) 
7 C 2 (V)+4C 2 (?A) -1 

Let us note that for any function G(^>) there are two particular solutions of 
Eq. (|61[). namely, 



Ci (-0) = « sin /??/>, C 2 (V>) = a cos /3ip : (69) 

where 

a = ±i /? = ±2. (70) 

Really, it is easy to check that the first particular solution C\ (ip) is the solution 
of Eq. (|6ip by direct calculations of the following terms: 

CiM+4CiM = (4-/3 2 )asin/#, (71) 

Cl{i]))+4el(i>)-l = (4 - p 2 ) a 2 sin 2 ^ + (o? F ~ l) • (72) 
Therefore, if we suppose that a and /? fulfil the following conditions: 

a 2 /? 2 -1 = 0, 4-/3 2 =0, (73) 

then for any function G(ip) the left- and right-hand sides of Eq. (j6"T1) are equal 
to zero separately. The same is true for the second particular solution C2(V0- 
If we differentiate Eq. ([77| with respect to ip, then we obtain that 

2Cty0(c'ty0+4C(7/0) =0. (74) 

If C(V>) = 0, then the function C(ip) is constant and from Eq. (|6"8")l we obtain 

GW= 4l^T W_W ' (75) 
If we take a non-constant (C(-0) ^ 0) solution 

C(V) = ±isin2(V-Vo) (76) 
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of the following differential equation: 

c 2 (V0 + 4C 2 = o, 



(77) 



then Eq. (j6Tj) is fulfilled identically for any function G(ip). 

The variables K and y can be expressed through C and its first derivative C 
as follows: 

2 



ctg^ 



C 



4C 



4C 2 



l+C 



1-AC 2 - C 2 



4C 2 



1 -C 



Really, using Eqs. (j3"9")l . (|52p . (I5"51) . and (|67l) we can calculate that 
Ksiny „ S 1 (l - AC 2 -C 2 ^ 



l+K 2 

Kcosy 
l + K 2 



= J = 



1 - 25 
C (1 + 2 J) = 



C 



i + 4C 2 +e 2 / 

2C 



K 2 = ^ '—. (78) 

(79) 
(80) 

K 2 + l ~^<-*" i_25 1 + AC 2 +C 2 ' (81) 

Solving Eq. (j8"Tj) with respect to 72. 2 provides us with the second of Eqs. (|78l) . 

Therefore, the complex admittance Y — Kexp(iy) can be expressed through 
C and its first derivative C as follows: 



1-25 l + AC 2 + C 2 



from where the first of Eqs. ([75)1 follows directly. 
Additionally, using Eq. (|62p we obtain that 

7? 2 — 1 . r 



2C 



Y = (l+K 2 



Kcosy Ksiny 
+ l 



l + K 2 



l + K 2 



AC + i 1 - 4C 2 - C : 



AC 2 



1-C 



(82) 



If the functions Q(tp) and C(ip) (or only one of them) are given, then the 
following expressions for X(ip) and the complex-valued wave function W (ip) can 
be written: 



W(ip) 



1 - G«W) 



di/> 



wq exp 



Wq exp 



i - G(ip)c(ip) ) y(^>)# 

1 -GC 



(83) 



4C + z 1 - 4C 2 - C' 



AC 2 



1-C 



dip 



(84) 



For the partial solutions of Eq. (jSTj). which are given by Eq. (|75|) . the ad- 
mittance Y(ip) is purely real, i.e., 



ctg (-0 - Vo) , 
-tg (tp - ip ) , 



(85) 
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therefore, for any given function G(tp) (equivalently Q(ip)) we obtain the fol- 
lowing expressions for X(ip) and the complex-valued wave function W(tp): 



W{ip) 



1 =F G{ip) sin (ip - ip ) cos (ip - ip ) 



dtp 



w 



V 7 ! - Z 2 W 



exp 



G(ip)Z 2 (iP)dip 



where 



(86) 
(87) 

(88) 
(89) 



where C(ip) is given by Eq. (|76[) , and for the complex admittance we obtain Eq. 
(|85p by direct calculation. Then Eq. (j59")l can be rewritten as follows: 



" + ": cos(^-^o), 
sin(^-Vo)- 



Really, to obtain Eq. (|87|l we need to calculate the following integral: 

\ - G(ip)C(ip)) Y 



Y(ip)dip- J G(tp)C(ip)Y(ip)dip = In Vl - ^ 2 W0 - / G(ip)Z 2 (%P)d%p, (90) 
where 

' + ; : sin ;^-^i' of) 

— : cos (t/> — t/> ) • 



Finally, we see that Eq. (JSTj) directly follows from Eqs. (JHOJ) and (|i?Tj) . 

For a constant function C(V0 and G('0) given by Eq. (1751) . we obtain the 
following expressions for Q(ip), the complex- valued admittance F(V0, -XXVOi 
and the complex- valued wave function W(ip): 



x{i>) 

W(ip) 



go exp 



8C 



4C 2 - 1 
AC + i (1 - 4C 2 

1 + 4C 2 
1 + 4C 2 1 
8C ' 
4C 

wo exp 



+ i \{ip - ip ) 



I -AC 2 

Eqs. (fM]) and (|95p directly follows from the fact that the expression 
and complex- valued admittance Y(ip), given by Eq. (j93| . are constants. 



(92) 

(93) 
(94) 
(95) 

(96) 
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5 Special solutions 



Though it is hardly possible to find the exact solution of Eq. (16TT) in a general 
case, the above analysis clarifies the nature of quasi-periodic Bloch waves in the 
transmission band and allows one to construct a wide class of special analytical 
solutions. A continual set of integrable wave equations can be obtained if we 
choose 

din M(C) _ 1 dM(C) 



GMC)] 



(97) 



dC M(C) dC 

where M(C) is an arbitrary real-valued function. In this case Eq. (f5T|) has an 
energy integral 



C 



?2 



1 - AC 1 + M(C) 



(98) 



and a periodic solution C(ip) = C(tp + t) given by the following expressions: 



t/j = ± 



dC 



y/1 - 4C 2 + M(C) 



dC 



(99) 



where the turning points C± are the roots of the radical. 

Let us notice that using Eq. (l97l) we can calculate the complete derivative 
of M [C(ip)] with respect to ip as follows: 



M(C) _ 1 dM{C) dC 
M(C) ~ M(C) dC dip 

Then rewriting Eq. (|6ip in the form 
2C(C + AC 



GC. 



GC 



In 



C 2 + AC' Z - 1 



C 2 + AC 2 - 1 dtp 
and integrating Eq. (jlOlj) we obtain that 

C 2 + 4C 2 - 1 



(100) 



lnM(C) (101) 



M{C) 



= const. 



(102) 



Due to the arbitrariness of the very function M(C), we can leave out the constant 
on the right-hand side of Eq. (|102[) . This leads us directly to Eq. (|98]) . 
Eq. (|61j) can be written in the following form: 



C = f[C,C,1> 



(103) 



where the direct dependence on -0 is realized only through the function G(jp). 
Let us suppose that in some way we have rewritten it as a function of C, i.e., 
G(C) = G[ip(C)\. Then in Eq. ([1031) tne direct dependence on ip is missing, 
and therefore, we can take C as an independent variable. Then we obtain that 
C = y(C), C = y(C)y'(C), and Eq. (JTD3]) reads 



2y(C)y'(C) - G(C)y 2 (C) = G(C) (AC 2 - l) - 8C, 



(104) 
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where / denotes the derivative with respect to C. Let us note that 



\M{C) J M{C) 



„ / M ( C ) 2 

yy M(C) y 



(105) 



and we see that to integrate Eq. (|104l) it is enough to suppose that the connection 
between the functions G(C) and M{C) is given by Eq. (|97p . Then we obtain the 
following first-order differential equation: 



C 2 = y 2 = M{C) 



AC 2 



M 2 {C) 



dM{C) 



8C 
M(C) 



dC 



(106) 



If we compare Eq. (|106[) with Eq. we obtain the compatibility condition 



4C 2 -1 
M(C) 



sc 



M 2 (C) 

If we suppose that the function M(C) is constant, i.e 
M(C) = c, c> -1, c^O, 
then G = and we obtain that 



M(C) 



dC. 



i>(C) = ± 



l 

c-4C 2 ~ 2 



2C 



y = 



(107) 



(108) 



(109) 



The integral in Eq. (I109P gives us the arcsin function, therefore, after the in- 
version we finally obtain that the auxiliary function C(ip) fulfilling Eq. (|6"Tj) has 
the following form: 



C(V0=±^^sin2(V>-^o): 



'00 = ^(O). 



(110) 



If we suppose that 

M(C) = c + 8eC, 
where c and e are constants, then 
dC 



c> -l-ie A 



V(C) = ± 



l 

= = ±- 

c + 8eC - 4C 2 2 



c^0, (111) 
2 (C - e) 



2/ = 



4e 2 ' 



(H2) 

After the integration and inversion of Eq. (|1 12[) with respect to the auxiliary 
function C(tp), we obtain that 



If we suppose that 



Vl + c + 4e 2 
e ± sin 2 \ip — ipo) . 



Inl 



M(C) = c + 8eC - d A C 



c > -1 



16e 2 
d 2 + 



c^0, 



(113) 



(114) 



14 



where c, e, and d are constants, then 



dC 



V = 



y/l + c + %eC - (d 2 + 4)C 2 
(d 2 + 4) C - 4e 



1 



dy 



V^V 2 



+ c) {d 2 + 4) + 16e 2 
After the integration and inversion of Eq. (|115[) we obtain that 

C(V) = ^— ^ {4e ± v/(l + c) (d 2 + 4) + 16e 2 sin [ d 2 + 4 - <0 O )] } 



(115) 
(116) 



If we suppose that 



(117) 



16e 2 



M(C) = c + 8eC+ (k 2 + 4)C 2 , c > -1 + -— , k > 0, c ^ 0, (118) 



where c, e, and A; are constants, then 



± 



dC 



±- 



1 



VI + c + 8eC + k 2 C 2 ~ 2 
k 2 C + 4e 



+ c)k 2 - 16e 2 



(119) 
(120) 



The integral in Eq. (|119l) gives us the hyperbolic arcsh function, therefore, after 
the inversion we finally obtain that 

C(i/j) = — |-4e ± V(l + c)fc 2 - 16e 2 sh fc - ^o)} ■ (121) 

Let us also consider an instructive example of modulated waves in a periodic 
dielectric medium, determined by the following potential: 



M(C) = (4a 2 - 1) + b 2 C 4 7 ab < 1, 4a 2 7^ 1, 



(122) 



where a and 6 are some constants (a, 6 > 0). Then the auxiliary function C(rp) 
is expressed through the Jacobi elliptic functions of the first kind [7] : 



C(V0 - ±ay/l+p 2 



2 - V>o) 



P= r ■ (123) 

ab 



Really, Inserting Eq. (|122l) into the first of Eqs. (IM1) we obtain that 



^ = , ± 1 

V4a 2 - 4C 2 + 6 2 C 4 h 



dC 



(C 2 - C 2 ) (C 2 - C 2 ) 



(124) 
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where the roots of the radical are given as follows: 

2 

&2 



Cl 



l±y/l 



b 2 



(125) 



If we take that C+ > C_ > C > (the roots C± are real for ab < 1; additionally 
the condition C+ > C_ imposes a6 7^ 1), then the integral in Eq. (|124j) can be 
expressed through the inverse Jacobi elliptic functions of the first kind: 



V(C) = ± T^sn- 



cT' el 



Let us denote 



Vl - v 7 ! - ^ 5 1 - Vl - a 2 6 2 



a 6 



Then inverting Eq. (|126[) and taking into account that sn(— a;) 
that 



V2 



a 2 & 2 = 

a 



„ „ 2a 



(126) 

(127) 
-sn x and 
(128) 



we finally obtain Eq. ([123]) . 

For any given function M[C(t/j)] the complex-valued wave function W(if)) 
can be rewritten as follows: 



W = wq exp 



± 



AC-iM 



(M - CM') dC 



AC 2 



(IT Vl-4C 2 + M) MVT^C 5 



M 



(129) 



Really, on the period we have to take the sign " +" while integrating from C_ to 
C+ and "— " while integrating back from C+ to C_ in the integral in Eq. (|129|) : 



AC - iM 



c- AC 2 

C- 



(1 - VI - AC 2 
AC - iM 



M) 



(M - CM') dC 

mvi-ac 2 Tm 

(M - CM') dC 



lc+ AC 2 + (1 + VI - 4C 2 + M) Wl - 4C 2 
Then we can rewrite Eq. (| 130[) as follows: 

AC - iM (M - CM') dC 



M 



(130) 



+ 



c_ 2 + M — 2Vl - AC 2 + M My/1 - AC 2 + M 
c + AC - iM (M - CM') dC 



2 + M + 2V1 - AC 2 + M My/1 - AC 2 + M 



(2 + M) 



4C 

/c_ (2 + M) 2 - 4 + 16C 2 - AM \ M 

And this leads us directly to Eq. (|132l) . 



(M - CM') dC 
y/l-AC 2 + M' 



(131) 
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In particular, for the complex increment we obtain that 



X + irj = In 



o f C+ 2 + Af f4C \ (M-CAI')dC 
J c _ M 2 + 16C 2 \m J Vl-4C 2 + M' 1 J 

Let us take an arbitrary even function M(C), then the function 

Gty) = I G(i>)di> = ± [ dMi ~ C) (133) 

i J M(C)^l-AC 2 + M(C) 

will be periodic. Moreover, for any even function M(C) we have x = Oj which 
means that |T / F( - 0)| is periodic, while the phase advance per period r, i.e., 

, = 4 / C+ 2 + M ( 13 4) 

determines the modulation period T — (2tt/t]) t of the quasi-periodic solution 
W(tp) predicted by the Floquets theory. 

For any even function M(C) the real part \ of the complex increment given 
by Eq. ([T32D, i.e., 

^ 2 + Af f CM'] 4CdC 



Af 2 + 16C 2 [ M J Vl - 4C 2 + M 



is an integral from C_ to C+ of the expression which is odd. Therefore, if we 
take the symmetric range of integration, then the integral of the odd expression 
vanishes and x = 0. The purely imaginary part r\ of the complex increment is 
an integral of the expression which is even, hence, taking the symmetric range 
of integration gives us exactly Eq. (|134l) . 
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